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ON  THE  DETERMINATION  OF 
HORIZONTAL  FORCES  A FLOATING  ICE  PLATE 
EXERTS  ON  A STRUCTURE 


Arnold  D.  Kerr 


INTRODUCTION  AND  STATEMENT  OF  PROBLEM 

When  a river,  lake,  or  part  of  a sea  freezes  over,  the 
ice  cover  formed  is  often  in  contact  with  a variety  of 
structures,  such  as  bridge  piers,  dams,  waterpower 
stations,  and  offshore  drilling  platforms.  Subsequently, 
any  change  in  the  position  of  the  ice  cover  (by  changing 
the  water  level,  by  moving  the  water  base  underneath 
it  or  the  air  above  it,  or  by  changing  the  air  tempera- 
ture) creates  forces  between  the  structures  and  the  ice 
cover  formed. 

Thus,  for  a rational  design  of  a structure  that  comes 
in  contact  with  a floating  ice  cover,  the  engineer  needs 
to  know  the  forces  the  ice  will  exert  on  this  structure, 
especially  the  largest  possible  force. 

Early  descriptions  of  problems  of  this  type  were 
given  by  Barnes6  and  by  Komarovskii.16  The  effort 
and  progress  made  in  this  field  during  the  following 
few  decades  were  very  modest,  as  evidenced  by  the  re- 
sults presented  at  the  American  Society  of  Civil 
Engineers  Symposium  on  Ice  Pressure  Against  Dams 
in  1954. 

During  the  past  decade,  the  problem  of  ice  forces 
has  gained  in  importance,  partly  because  of  the 


discovery  of  oil  in  the  northern  regions  and  the  result- 
ing increased  activity  in  this  area.  The  intensified  re- 
search effort  that  has  followed  is  evidenced  by  the 
recent  increase  of  conferences  on  this  subject,  such  as 
the  Conference  on  Ice  Pressures  Against  Structures8 
held  in  1966,  the  Ice  Seminar1 1 held  in  1968,  the  Inter- 
national Conferences  on  Port  and  Ocean  Engineering 
under  Arctic  Conditions13  held  in  1971,  1973  and  1975, 
and  the  IAHR/AIRH  Symposia  on  Ice  and  Its  Action 
on  Hydraulic  Structures12  held  in  1970, 1972,  and 
1975.  In  spite  of  all  these  efforts,  to  date  no  reliable 
methods  have  been  developed  for  predicting  the  ice 
forces  an  ice  cover  may  exert  on  a structure. 

The  analyses  of  these  ice  forces  are  usually  grouped 
into  several  categories,  for  example:  1 ) analyses  of 
forces  caused  by  a rise  or  drop  in  the  water  level, 

2)  analyses  of  forces  caused  by  a structure  constraining 
the  movement  of  the  ice  cover  in  the  horizontal  plane, 
and  3)  analyses  of  forces  caused  by  the  impact  of  a 
moving  ice  flow  and  a structure.  This  report  discusses 
problems  of  the  second  category.  Examples  of  such 
problems  are  shown  in  Figure  1 . 

Since  ice  is  a viscoelastic  material,  the  forces  an  ice 
cover  exerts  on  a structure  are  time  dependent.  The 


a.  ice  cover  edge  pressure  b.  Ice  covet  pressing  against  an  isolated 

against  rigid  wall.  structure. 

Figure  1.  Structures  which  constrain  the  movement  of  the  ice  cover. 


Figure  2.  Ice  cover  constrained  by  two  parallel  rigid  walls,  p = ice  pressure 
against  wall  (per  unit  length  of  wall)  and  h = thickness  of  ice  cover. 


study  of  these  forces,  caused  by  constrained  thermal 
expansions,  was  initiated  by  Royen  in  1922.  The 
analytical  and  test  results  of  Royen,  as  well  as  the 
more  recent  results  of  a number  of  other  investigators, 
were  reviewed  by  Korzhavin17  (chapter  XIV)  in  1962, 
by  Drouin7  in  1966,  and  by  Michel18  in  1970.  Other 
attempts  to  solve  this  problem  were  presented  by 
Nevel19  in  1968  and  by  Panfilov21  in  1965. 

The  conclusion  of  the  survey  by  Drouin7  (still  valid 
presently)  was  that  the  determination  of  ice  forces  as 
a function  of  time  is  not  yet  solved  satisfactorily.  The 
relevant  results  of  references  19  and  21  are  not  con- 
clusive either.  A major  difficulty  in  solving  this  prob- 
lem is  the  absence  of  appropriate  constitutive  equations 
for  the  ice  cover.  It  should  also  be  noted  that,  even  if 
expressions  which  determine  the  ice  forces  were  avail- 
able, the  magnitude  of  these  forces  may  sometimes  be 
limited  by  a possible  failure  (crushing  or  buckling)  of 
the  ice  cover.  This  aspect  will  be  discussed  later. 

For  the  design  of  a structure  that  is  in  contact  with 
an  ice  cover,  it  is  often  sufficient  to  know  the  largest 
force  the  ice  cover  may  exert  on  the  structure.  This 
approach  simplifies  the  necessary  analyses. 

To  demonstrate  the  essential  aspects  of  this  ice  force 
problem,  consider,  as  an  example,  an  ice  cover  con- 
strained by  two  parallel  "rigid”  walls,  as  shown  in 
Figure  2. 

Assume  that  (at  a time  / = 0)  the  temperature  of 
the  ice  plate  is  instantaneously  raised  by  a constant 
value,  AT°.  Then,  the  corresponding  axial  compression 
force  p0  (per  unit  length  of  wall)  is  approximately 

p0=EhaAT°.  (1) 

Equation  1 is  based  on  the  assumption  that  at  the 
instant  of  the  temperature  rise  (t  = 0)  the  ice  obeys 
Hooke’s  law,  that  Young’s  modulus  E is  a value  aver- 
aged across  the  plate  thickness,  and  that  the  coefficient 
of  linear  expansion  a is  also  an  average  value  which 
doesn’t  vary  noticeably  during  the  temperature  in- 
crease AT°. 


Because  of  the  viscoelastic  response  of  the  ice  cover, 
the  pressure  p,  for  t > 0,  is  smaller  than  the  p0  value 
given  in  eq  1 : 

P(t)  < p0.  (2)  > 

A similar  argument  leads  to  the  conclusion  that 
eq  2 is  also  valid  for  the  case  when  the  temperature  rise 
is  not  instantaneous,  but  proceeds  gradually,  as  is  the 
case  in  the  field.  Thus,  eq  1 yields  the  largest  possible 
axial  compression  force  due  to  AT0.  This  suggests 
that,  for  many  engineering  problems,  it  may  be  suffi- 
cient to  determine  the  forces  acting  on  a structure  by 
assuming  that  the  ice  cover  responds  elastically. 

At  this  point,  it  should  be  noted  that  the  force  esti- 
mate given  by  eq  1 , or  the  more  accurate  expression 
which  takes  into  consideration  the  viscoelastic  response 
of  the  ice  cover,  is  valid  as  long  as  the  ice  cover  does  not 
fail.  Thus,  beyond  a certain  compression  force,  eq  1,  or 
its  viscoelastic  equivalent,  yields  too  high  values  for  the 
ice  force  and  it  has  to  be  bounded  by  a failure  criterion. 
This  criterion  is  obtained  by  noting  that  an  ice  force 
cannot  be  larger  than  the  force  capable  of  breaking  up 
the  ice  cover.  Thus,  the  needed  bound  for  in-plane  com- 
pression forces  is  established  by  the  determination  of 
the  intensity  of  the  force  at  which  the  ice  cover  fails. 
Note  that  this  criterion  is  valid  for  the  determination  of 
the  largest  ice  forces  that  may  be  due  to  a change  in 
temperature  or  due  to  any  other  reason,  such  as  the 
movement  of  the  water  in  the  horizontal  direction. 

Many  published  analyses  are  based  on  this  notion. 

In  the  majority  of  these  publications,  it  is  assumed  that 
this  failure  load  is  directly  related  to  a "crushing 
strength”  of  the  ice;  that  is,  the  mechanism  of  failure  is 
due  to  crushing  or  splitting  of  the  floating  ice  plate  in 
the  immediate  vicinity  of  the  structure.  A review  of  the 
relevant  literature  up  to  1962  is  given  by  Korzhavin.17 
A more  recent  use  of  this  approach  is  presented  by 
Afanasiev,1  2 Assur,  5 Michel,18  Schwarz,25  and  Shadrin 
and  Panfilov,26 
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In  laboratory  and  field  tests  it  was  observed,  how- 
ever, that,  for  relatively  thin  ice  plates  (compared  with 
the  width  of  contact  of  plate  and  structure),  the  ice 
cover  failed  by  buckling  in  the  vicinity  of  the  contact 
area.3  1 5 20  24 

The  published  analyses  of  ice  forces  on  isolated 
structures,  which  are  based  on  the  buckling  failure 
mechanism,  are  not  conclusive,  however.  The  purpose 
of  this  report  is  first  to  review  these  few  attempts  and 
then  to  develop  an  improved  method  of  analysis  for 
the  determination  of  the  largest  force  that  a relatively 
thin  ice  cover  may  exert  on  an  isolated  structure.  The 
proposed  analysis  is  based  on  the  buckling  mechanism 
of  the  floating  ice  plate. 


REVIEW  OF  RELEVANT  ANALYSES 
AND  TESTS 

An  early  attempt  to  determine  analytically  the 
buckling  force  for  an  ice  plate  was  presented  by 
Albenga4  in  1921.  However,  because  he  neglected  the 
effect  of  the  liquid  base,  he  obtained  results  that  do 
not  apply  to  floating  ice  covers. 

On  the  basis  of  observations  in  the  field,  in  1 947 
Rose24  (p.  883)  stated:  "...  that  arching  and  buckling 
of  an  ice  sheet  does  occur  but  that  there  is  a limiting 
thickness  of  about  1 ft  beyond  which  buckling  will 
rarely  occur.”  He  also  claimed  (p.  880):  ”...  for  ice 
sheets  less  than  1 ft  thick,  the  critical  buckling  load 
is  the  limiting  factor  governing  the  ice  pressure  on  a 
dam,  whereas  for  very  thick  ice  sheets  (about  2 ft  or 
more),  the  crushing  strength  pressure”  is  the  limiting 
factor.  No  analysis  was  given  for  the  determination  of 
the  buckling  load. 

In  1962,  Shadrin  and  Panfilov26  presented  an 
analysis  for  the  semi-infinite  ice  cover  whose  straight 
edge  is  pressing  against  a rigid  wall,  as  shown  in 
Figure  la  and  Figure  2.  The  authors  assumed  that  the 
type  of  ice  cover  failure  depends  on  the  thickness  of 
the  ice  cover,  namely,  buckling  for  thin  covers  and 
local  crushing  at  the  contact  area  for  thick  covers. 

This  assumption  is  in  agreement  with  the  observations 
made  by  Rose,24  as  described  above. 

Shadrin  and  Panfilov26  expressed  the  failure  load 
(per  unit  length)  as 

Pi  = °ct h (3) 

where  ocf  is  a compression  strength  and  h is  the  thick- 
ness of  the  ice  cover. 

For  relatively  thin  covers,  the  buckling  loadpb  was 
determined  from  the  differential  equation 


D + p + yw  = 0 (4) 

dx 4 dx 2 

and  the  boundary  conditions 
w(0)  = 0 

■2  <5> 
£2-  =0 

dx 2 o 2 

where  w(x)  is  the  vertical  deflection  at  x and  at  any 
value  -»<  y < °°,  D = £7j3/[12(1  - r2)]  is  the  flexural 
rigidity  of  the  cover,  and  y is  the  specific  weight  of  the 
liquid  base.  The  solution  of  the  above  formulation 
yields 

Pb  = 2\/yD  . (6) 

The  above  expression  forpb  is  also  valid  when  the 
boundary  at*  = 0 is  clamped,9  as  may  often  be  the 
case  in  the  field. 

Next,  eq  3 and  6 are  compared  by  noting  that, 
according  to  tests,  if  the  corresponding 

Pb<Pi 

the  ice  cover  fails  by  buckling  and  if 


Pi  <Pb 

the  ice  cover  fails  by  crushing.  Substituting  eq  3 and  6 
in  the  first  inequality,  it  follows  that  the  ice  cover  will 
fail  by  buckling  when 


h < 


3(1  -v2)  0% 

yE 


(7) 


As  an  example,  if  v = 0.34,  y = 0.001  kg/cm3,  E = 
40,000  kg/cm2,  and  ocf  = 20  kg/cm2,  then  according 
to  eq  7 the  ice  cover  will  buckle  when 


/k  27  cm  s 10.5  in. 


This  result  agrees  reasonably  well  with  the  field  observa- 
tions reported  by  Rose.24 

According  to  the  above  analysis,  when  inequality  (7) 
is  satisfied,  the  largest  force  that  the  ice  cover  may  exert 
on  the  long  rigid  structure  is  the  buckling  loadpb  = 
2\/yD\  otherwise,  the  largest  force  isp(  = ocf/t. 

Next,  the  calculation  of  the  largest  force  an  ice  cover 
exerts  on  a ship’s  hull,  presented  by  Kheishinls  in  1961 
is  reviewed;  this  is  a problem  closely  related  to  the  deter- 
mination of  ice  forces  on  isolated  structures.  Kheishin 
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Figure  3.  Analytical  mode / used  by  Kheishin. 


reduced  this  problem  to  the  determination  of  the  pres- 
sure p = constant  at  which  a floating  elastic  semi- 
infinite plate,  shown  in  Figure  3,  will  fail.  He  based 
his  analysis  on  field  observations  that  ice  cover  failures 
often  occur  in  two  stages:  the  formation  of  cracks 
which  emanate  from  points  1 and  2 (Fig.  3),  followed 
by  the  buckling  of  the  inner  part  of  the  ice  cover, 
which  is  formed  by  these  cracks. 

For  his  analysis,  Kheishin1  5 assumed  that  the 
cracks  are  initiated  by  the  largest  shear  stress  along 
the  edge 

KyLax  = axy  (°,  ± */2)  = pM.  (8) 

Choosing  for  the  shear  strength  of  the  (sea)  ice  cover 
osf  = 30  tons/m2,  and  using  the  equation 

(°xy)max  - °st  (9) 

as  a crack  criterion,  he  obtained  from  eq  8 the  pressure 
intensity  that  causes  the  formation  of  the  cracks  as 


Pjf  s 94/7  tons/m 


(10) 


where  h is  in  meters. 

To  simplify  the  buckling  analysis,  Kheishin15  as- 
sumed that  the  two  cracks  which  emanate  from  the 
points  (0,  ± b/ 2)  are  parallel  to  the  x-axis  and  that  the 
formed  ice  strip  buckles  when 


P b 


= VyD . 

cyl  bend 


(ID 


Assuming  v = 0.36,  y = 1 ton/m3,  and  £ = 40,000 
tons/m2,  the  above  equation  reduces  to 

pb|  =62  h\/h  tons/m  (12) 

I cyl  bend 

where  h is  in  meters. 


* Note  that  in  the  previous  example  (Shadrin  and  Panfilov1* ) £ 
interval  is  0 < h < 0.23  m. 


Comparing  eq  10  and  eq  12,  Kheishin15  obtained 
that  for  0 < h < 2.2  m,*  which  is  the  range  of  practical 
interest,  the  inequality 

^Hylbend  (13) 

holds  (although,  according  to  Kheishin,15  p^  < pb, 
pb  being  the  buckling  load  of  the  uncracked  plate). 

Based  on  observations  in  the  field,  that  an  ice  cover 
cracks  prior  to  buckling,  Kheishin  concluded  that  the 
largest  ice  force  is  given  by  pit,  in  eq  10.  According 
to  the  above  analysis,  after  the  plate  cracks,  buckling 
of  the  inner  strip  takes  place  immediately,  since 

Pb|cylbend<psC  04) 

The  above  method  is  based  on  a number  of  question- 
able assumptions,  such  as  that  eq  9 is  a valid  criterion, 
that  the  cracks  are  parallel  and  normal  to  the  free  edge, 
etc.  Note  also  thatpb  given  in  eq  1 1 is  the  critical 
buckling  pressure  when  the  edge  is  free  to  rotate  and 
displace  vertically,9  which  is  usually  not  the  case  when 
the  ice  cover  is  in  contact  with  a fixed  vertical  structure 
or  with  the  hull  of  a ship.  Other  shortcomings  of  eq  1 1 
are  that  the  ice  strip  formed  by  two  cracks  that  emanate 
from  points  1 and  2 is  not  of  constant  width  and  that 
the  strip  is  constrained  along  the  entire  length  of  the 
formed  edges  by  the  surrounding  ice  cover.  Also, 
according  to  Nevel  ct  al.,20  who  conducted  many  closely 
related  tests,  the  formation  of  radial  cracks  is  not  followed 
immediately  by  ice  cover  failure.  Therefore,  the  results 
obtained  by  Kheishin,15  which  are  also  reproduced  in 
the  book  by  Popov  et  al.,23  should  be  used  with  caution. 

In  1971,  Afanasiev  et  al.3  published  results  of  ar 
extensive  laboratory  test  program  in  which  a number 
of  thin,  floating  ice  plates  (2.6  § h $ 3.4  cm)  were  sub- 
jected to  a horizontal  force  by  pressing  a circular  or  a 
rectangular  pier-model  against  the  free  edge  of  the  plates. 
According  to  their  results,  when  b/h  > 4 and  the  ice 
cover  was  not  made  more  rigid  by  the  proximity  of  the 
basin  walls  or  the  presence  of  simulated  ridges,  the  ice 
plate  failed  by  buckling.  During  these  tests,  it  was  ob- 
served (ref.  3,  p.  64)  that,  at  a certain  intensity  of  the 
load,  radial  cracks  formed  in  the  plate.  At  a further 
increase  of  the  load,  the  ice  plate  lifted  along  a nearly 
circular  curve  at  a distance  of  about  1 4 x/z  from  the 
pier.  This  was  followed  by  a breakdown  of  the  ice  field. 

Afanasiev  et  al.3  then  presented  an  analysis  for  the 
semi-infinite  ice  plate  pressing  uniformly  against  a 
rigid  wall,  which  is  similar  to  the  one  given  by  Shadrin 
and  Panfilov26  discussed  previously.  The  results  they 
obtained  arc  also  similar:  for  fresh  water  ice  the  ice 
cover  buckles  when  h < 30  cm,  and  for  salt  water  ice 
r 400,000  tons/m2  was  used.  For  this  E -value  the  corresponding 
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a.  Rectangular  pile. 


Figure  4.  Observed  crack  patterns  for  different  shapes  of  model  pier. 


when  h < 20  cm.  Thus,  when  these  inequalities  are 
satisfied,  the  buckling  load  is  the  largest  force  the  ice 
cover  can  exert  against  a rigid  wall. 

In  1972,  Nevel  et  al.20  presented  results  of  exten- 
sive ice  force  tests  conducted  at  CRREL.  The  test 
set-up  was  similar  to  the  one  described  by  Afanasiev 
et  al.,3  except  that  the  basin  used  (6.3x6 .3  m)  and  the 
thickness  of  the  plates  tested  (6.5  S h rS  22.4  cm) 
were  larger.  Their  description  of  the  test  results  is 
also  fuller.  The  results  of  tests  with  rectangular  piers 
showed  that  for  an  average  ice  temperature  of  -4.3°C 
and  h =*  12  cm,  for  b/h  > 2 the  ice  plate  failed  by 
buckling.  For  circular  piers,  instead  of  the  6-value, 
they  used  the  pier  diameter  d.  From  the  test  results 
presented,  it  follows  that  for  cold  ice  (-4.0°C)  the  ice 
plate  buckled  for  d/h  > 4.85,  whereas  for  warmer  ice 
(-2.7°C)  the  ice  plate  buckled  for  larger  ratios  of  d/h: 
d/h  > 8.2.  In  some  of  these  tests,  it  was  observed  that, 
as  the  load  increased,  cracks  occurred,  as  indicated  in 
Figure  4.  For  rectangular  piles  (Fig.  4a),  the  cracks 
formed  a funnel-shaped  plate  which  buckled  in  the 
wider  part  of  the  plate. 

According  to  the  above  test  results,3  20  for 
rectangular  piers,  a thin  ice  plate  may  fail  in  buckling 
for  b/h  > 2 and  for  circular  piers  for  d/h  > 5,  ranges 
of  practical  interest.  Thus,  for  a variety  of  isolated 
structures  and  relatively  thin  ice  covers,  the  largest  ice 
forces  appear  to  be  related  to  the  buckling  load  of  the 
ice  cover  and  not  to  the  local  "crushing  strength.” 

The  results  of  the  reviewed  analyses  presented  by 
Shadrin  and  Panfilov26  and  by  Afanasiev  et  al.3  for  the 
determination  of  the  largest  force  an  ice  cover  may 
exert  on  a long  vertical  wall  (for  example,  a long  dam) 
are  in  agreement  with  observations  in  the  field,  regard- 
ing occurrence  of  buckling  failures.  Therefore,  it  is 
reasonable  to  expect  that  if  the  material  parameters 


are  chosen  properly,  the  largest  ice  force  for  problems 
of  this  type  can  be  estimated  using  the  analyses  dis- 
cussed: eq  4 to  7. 

However,  for  the  determination  of  the  largest  ice 
forces  on  isolated  structures,  based  on  the  buckling 
failure  mode,  only  the  analysis  of  Kheishin1  s seems 
to  be  available.  As  shown  above,  his  analysis  is  based 
on  questionable  assumptions.  According  to  Kheishin, 
the  failure  load  is  the  load  at  which  the  plate  cracks, 
since  buckling  takes  place  immediately  after  the  cracks 
form;  but  this  is  not  substantiated  bv  the  observations 
described  by  Afanasiev  et  al.3  and  Nevel  et  al.20 

These  findings  suggest  the  need  for  an  improved 
analysis  for  the  determination  of  the  largest  forces  an 
ice  cover  exerts  on  isolated  structures,  based  on  the 
buckling  failure  mechanism.  Such  an  analysis  is  pre- 
sented in  the  following  section. 

DETERMINATION  OF  THE  LARGEST 
ICE  FORCE  ON  AN  ISOLATED  STRUCTURE 

Preliminary  remarks 

According  to  the  laboratory  tests  by  Nevel  et  al.,20 
when  a pier  with  a flat  surface  presses  against  the  free 
edge  of  a floating  ice  plate  (or  vice  versa),  at  a certain 
intensity  of  the  in-plane  load  radial  cracks  emanate 
from  the  loaded  region,  as  indicated  in  Figure  4a; 
a further  increase  of  the  load  leads  to  more  radial 
cracks  and  finally  to  the  buckling  of  the  cracked  region. 
A similar  observation  was  also  made  by  Afanasiev  et 
al.3  These  observations  suggest  that  the  largest  force  a 
relatively  thin  ice  cover  could  exert  on  an  isolated 
structure  may  be  related  to  the  buckling  load  of  float- 
ing wedge-shaped  plates. 
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The  buckling  analysis  of  a floating  wedge 

Consider  a floating  truncated  wedge,  as  shown  in 
Figure  5.  It  is  assumed  that  the  buckling  load  Pb  may 
be  determined,  approximately,  from  the  eigenvalue 
problem,  consisting  of  the  differential  equation 

[£/(x)w"|  + Pw"  + k(x)w  = 0 0 < x < °°  (15a) 

the  boundary  conditions  at  the  contact  area 


w(0)  = 0 
w"(0)  =0 


(16a) 


or 


w(  0)  = 0 


w'(0)  =0 


(16b) 


and  the  regularity  conditions  at  x = °°. 

In  the  above  formulation,  w(x ) is  a vertical  pertur- 
bation of  the  wedge  axis  at  x.  Furthermore, 


M*)  = y(t>0  + x<p*) 
El(x)  = D(b0  + x0*) 


(15b) 


where 


0 = gftl 

12(1  -i>2) 

<p*  = 2 tg-£ 


(15c) 


and  y is  the  specific  weight  of  the  liquid  base.  The 
term  (1  -f2)  was  inserted  in  the  denominator  of  £/(x) 
to  account  for  plate  action,  as  done  by  Papkovich22 
in  connection  with  a related  problem,  recently  dis- 
cussed by  Kerr14  (p.  15). 

Substituting  eq  1 5b  into  eq  1 5a,  results  in 

[D(b0+xO*)w"\  + Pw" 


+ y(b0  + x<p*)w  = 0 0 < x < °°  (15d) 

a fourth  order  equation  with  variable  coefficients. 

For  <t>*  = 0 (i.e.,  for  <f>  = 0),  the  above  equation  reduces 
to  the  equation  with  constant  coefficients  used  by 
Kheishin,15  which  is  the  same  as  eq  4. 

Because  of  the  nature  of  the  variable  coefficients, 
no  exact  close  form  solution  is  available.  A numerical 
solution  can  be  easily  obtained  using,  for  example,  a 


F igure  5.  Floating  truncated  ice  wedge. 


finite  difference  scheme,  because  the  buckling  mode 
is  restricted  to  the  vicinity  of  the  contact  region  at 
x = 0.  However,  the  solution  obtained  will  be  valid 
only  for  the  assumed  set  of  parameters. 

For  the  proposed  method  to  determine  £max,  an 
analytical  expression  that  shows  the  dependence  of  the 
buckling  load  Pb  on  the  various  geometrical  and 
material  parameters  of  the  floating  ice  wedge  and  the 
width  of  the  structure  is  more  suitable.  Noting  that 
the  representation  of  the  buckling  of  a cracked  ice 
cover  by  the  buckling  of  a floating  wedge  is  at  best  an 
approximation,  it  appears  sufficient  to  solve  the  wedge 
formulation  approximately.  In  the  following,  such  an 
analytical  solution  is  obtained,  at  first,  using  the 
Stodola-Viancllo  ,nc'>  >d.10 

For  this  purpose,  eq  15d  is  rewritten  as 

w,v  + 4K4w  = -|lj(xw1")  +2-xwi|-X2w1”  (17a) 


where 


k4  = EX— 
4 D 


X2 


_ JP_ 
Dbn 


Assuming  for  Wj  the  expression 
w,  (x)  = A , c~KX  sin  kx 


(17b) 


(18) 


and  substituting  it  into  eq  17a,  a nonhomogeneous 
ordinary  differential  equation  with  constant  coefficients 
for  w(x)  is  obtained.  Its  solution,  which  satisfies  the 
boundary  conditions  of  a simply  supported  edge,  eq  16a, 
is 


w(x)  - A , e*“*  | A sin  (xx)  + (xA  +8)x  sin  (xx) 

+ fix  cos(xx)]  (19a) 
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where 


A = 


2Kb0 


. P 
&Dxb0 


> 


(1%) 


The  first  eigenvalue  is  obtained  by  substituting 
eq  19a  into  the  condition 


as  an  approximate  shape  function,  substituting  it  into 
V2,  and  performing  the  necessary  integrations,  V2  be- 
comes 

V2  =-j^{(6Db0p4  -Pp2  +yb0/2)p 

+ 4>*(2Dp4  +7/2)| . (24) 

If  P is  chosen  a priori  as  the  value  which  corresponds 
to  the  compressed  semi-infinite  plate, 


J w\[x)dx  = J w](x)w(x)dx.  (20) 

0 0 

Performing  the  necessary  integrations  and  noting 
that  X2  = P/(Db0),  the  approximate  buckling  load  is 
obtained  as 

(/>b)sirnply=5-3OK(K6o+0*).  (21a) 

supp 


p-  k'  = \/yjD  - \/4k  , (25) 

it  follows  that  V2  = V2(A).  The  condition  6 V2  = 0 then 
reduces  to 


(26) 


The  above  equation  yields 


If,  instead  of  the  boundary  conditions  for  a simply 
supported  end,  the  ones  for  a clamped  end  at  x = 0 
are  used,  the  above  procedure  yields 


(/5b)clamped=8°K(2^o  + «*) 


(21b) 


Note  that  in  both  cases  the  respective/^  expression 
consists  of  the  sum  of  two  terms: 


{6Db0K4  + yb0/2  - PbK'2)n' 

+ <t>*  (2Dk'4  + y/2)  = 0.  (27) 

Thus 


(DK2b0) 

(Dk<P*) 

each  multiplied  by  a constant  coefficient.  Thus,  for 
the  obtained  Pb  expressions,  the  different  boundary 
conditions  at  x = 0 affect  only  these  coefficients. 

To  check  if  the  above  two  terms  are  representative 
of  the  (bifurcation)  buckling  solution  for  the  com- 
pressed wedge,  in  the  following  the  wedge  problem  is 
sdved  again  using  the  energy  method. 

With  the  notation  used  above,  the  second  variation 
of  the  total  potential  energy  for  the  wedge  problem 
shown  in  Figure  5 is 


V2=f(f  W"2  + J w2~t  w'2)dx  (22) 

where  El  and  k are  as  given  in  cq  15b. 

Assuming 

w(x)  = A e"®*  sin/3x  (23) 


and  noting  that  7/(1 2 Dk'4)  =1/12  and  k'  = \/4  k,  Pb 
reduces  to 


Pb  = 3.5Dk  [ 3.760k  + 0* | . (28b) 


If  P is  not  fixed  a priori,  then  V2  = V2  (A,  P).  The 
condition  6 V2  = 0 then  reduces  to 


bV2 

W 


= 0. 


(29) 


The  first  equation  in  (29)  is  identical  to  eq  26,  and  yields 
eq  28a  with  k'  = P: 


Pb=6Dp\bQp(\  + 


(30) 


The  second  equation  in  (29),  after  eliminating  Pb  by 
using  eq  30,  yields 
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Figure  6.  Critical  buckling  loads  for  float- 
ing wedge.  (Note:  Pb  = b0  2\fyD  ). 


P(4Db0p4  - b0yl 3)  + <t>*  (2D/34/3  - y/2)  = O.f  (31 ) 

This  is  the  equation  for  the  determination  of  (3  which 
appears  in  eq  30  as  an  unknown. 

Comparison  of  eq  28b  with  eq  21a  and  eq  21b  shows 
that,  except  for  the  constant  coefficients,  the  three 
equations  have  the  same  functional  form. 

In  order  to  compare  the  approximate  expressions 
obtained  for  the  buckling  load  Pb,  they  were  evaluated 
for  E = 30,000  kg/cm2,  v = 0.34,  y = 0.001  kg/cm3, 
b0  = 20  cm,  and  h = 10  cm.  The  results  are  shown  in 
Figure  6. 

Proposed  method  to  determine  Pmlx 

When  an  ice  cover  presses  against  a bridge  pier  or 
the  leg  of  a drilling  platform,  the  boundary  conditions 
at  the  contact  area  between  the  ice  and  the  structure 
fluctuate,  depending  on  a number  of  factors.  Also, 
the  value  <p  is  usually  not  a measurable  entity,  since 
according  to  test  observations,  at  the  onset  of  buckling, 
the  ice  cover  exhibits  many  radial  (and  some  nonradial) 
cracks. 

It  is  therefore  proposed  to  express  the  largest  force 
an  ice  cover  exerts  on  an  isolated  structure  (when 
buckling  dominates)  by  the  expression 

x=n'DK(Kb0+m')  (32) 

where,  according  to  eq  15c  and  eq  17b, 


f Note  that  eq  31  may  be  obtained  directly  from  the  condition 
»Pbl»0  = 0. 


£>  = _i*2 

12(1  -1>2) 


K = VyJ*D 


(33) 


and  n'  and  m'  are  parameters  to  be  determined  from 
field  tests. 

Since  the  effective  values  of  £ and  v depend  upon 
the  nature  of  the  precracked  plate,  the  term 
£/[12(1  -i'2)]  = E*  should  be  determined  from  the 
same  test  data.  Thus,  the  proposed  expression  for 
^max  *s 


Pmax  = n i/yE*W  W(-y/£*/t3)  b0  + m ] (34) 


where  y is  the  specific  weight  of  the  liquid  base,  h is  the 
ice  cover  thickness,  and  b0  is  the  width  of  pier.  The 
parameters  m,  n,  and  E*  are  to  be  determined  from 
appropriate  test  results. 

It  is  expected  that  m , n,  and  E*  will  vary  for  differ- 
ent categories  of  problems  and  different  ice  temperatures 
(in  particular  in  the  vicinity  of  the  melting  temperature). 
This  point,  as  well  as  the  validity  of  eq  34  to  represent 
/’max  (when  the  buckling  failure  mode  predominates), 
can  only  be  established  by  comparing  eq  34  with  the 
relevant  results  of  laboratory  and  field  tests. 
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In  accordance  with  letter  from  DAEN-RDC , DAEN-ASI  dated 
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Laboratory  Technical  Publications,  a facsimile  catalog 
card  in  Library  of  Congress  MARC  format  is  reproduced 
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On  the  determination  of  horizontal  forces  a floating  ice 
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